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Abstract 

Metastable, supersymmetry-breaking configurations can be created in flux geometries by placing 
antibranes in warped throats. Via gauge/gravity duality, such configurations should have an inter- 
pretation as supersymmetry-breaking states in the dual field theory. In this paper, we perturbatively 
determine the asymptotic supergravity solutions corresponding to D3-brane probes placed at the tip 
of the cascading warped deformed conifold geometry, which is dual to an SU(N+M) x SU(N) gauge 
theory. The backreaction of the antibranes has the effect of introducing imaginary anti-self-dual 
flux, squashing the compact part of the space and forcing the dilaton to run. Using the generaliza- 
tion of holographic renormalization to cascading geometries, we determine the expectation values 
of operators in the dual field theory in terms of the asymptotic values of the supergravity fields. 



1 Introduction 



There are a number of reasons to be interested in theories which break supersymmetry at an expo- 
nentially low scale. They may be of use as hidden sectors in supersymmetric models of weak-scale 
physics. The dynamics involved in the supersymmetry breaking is also often intricate, exhibiting 
novel phenomena in field theory. 

One relatively new tool to study such theories is gauge/gravity duality [lj. Type IIB string 
theory on spaces of the form AdSs x X 5 with X 5 a five- dimensional Sasaki-Einstein space is known 
to be dual to a class of four- dimensional superconformal field theories [2]. Simple examples with 
varying amounts of supersymmetry can be generated as orbifolds of S 5 [3] , or starting from conical 
hypersurface singularities like the conifold [4]. The dual field theories can be constructed by studying 
the gauge theories on N D3-branes at the appropriate singularities. 

Of more interest to us will be the backgrounds which are dual to confining M = 1 gauge theories 
in four dimensions. A canonical example is provided by the theory of M fractional D5-branes and 
iV D3-branes on the conifold, with M <C N '. This theory exhibits a "cascade" of Seiberg dualities 
[5], ending with chiral symmetry breaking at an exponentially small scale. The dual supergravity 
solution, found by Klebanov and Strassler [B] , is completely smooth; the "tip" of the AdS-like throat 
is rounded off by a complex deformation, whose magnitude is related to the scale of chiral symmetry 
breaking. 

It was proposed in [7] that one can make metastable supersymmetry breaking states starting 
from such a geometry, by adding a small number p <C M of D3-branes to the background. In 
the probe approximation, one can study the dynamics by using the DBI action on the antibranes. 
They are attracted to the smooth tip of the geometry by the background 5-form flux and warping 
of the metric. The D3-branes undergo the Myers effect |8J, but for sufficiently small p/M they 
cannot perturbatively decay to a supersymmetric vacuum. Instead, the leading decay channel 
is via quantum bubble nucleation, which is highly suppressed for suitable choices of parameters. 
Because the warp factor at the tip is exponentially small, these metastable states naturally break 
supersymmetry at an exponentially small scale, and thus should be usefully thought of as string 
duals of dynamical supersymmetry breaking. These supersymmetry-breaking vacua are still largely 
mysterious from the field theory point of view. 

For many purposes, it would be useful to have a supergravity solution which captures the 
long-range fields of the D3-brane perturbation to the Klebanov- Strassler geometry. Most notably, 
in gauge/gravity duality, one reads off the properties of states in a field theory by studying the 
(normalizable modes of) background supergravity fields. So if one wishes to compute the vacuum 
energy and the VEVs of the operators in this supersymmetry-breaking state of the Klebanov- 
Strassler field theory, one should find the supergravity solution in the presence of the D3-branes 
and read off these quantities from the behavior of various fields at infinity^ 

In this paper, we find the supergravity solutions characterizing these states, in perturbation 
theory in p/N and far from the tip of the throat. More precisely, we linearize the IIB equations 
of motion around the large r limit of the Klebanov- Strassler solution (first found by Klebanov and 
Tseytlin [IS]), and identify the modes which are activated by the presence of the p D3-branes. 
To make the computations more tractable, we study a "smeared" solution, where we imagine the 

*A distinct class of non-supersymmetric vacua has been explored in [9], building on the work of [lOj . The theory 
investigated in [9] involves a perturbation of the Klebanov-Strassler field theory Lagrangian by a non-supersymmetric 
gaugino mass term, which maps to a non-normalizable supergravity perturbation. In contrast, we are studying states 
of the Klebanov-Strassler theory itself, and only consider normalizable perturbations on the gravity side. Other work 
exploring similar issues in closely related theories appears in [TTJ [T21 [TJl Q3] . 
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D3-branes are spread uniformly over the internal dimensions at the tip. 

The result is simple to describe: in the supersymmetric theory, the geometry of the internal 
dimensions is the Sasaki-Einstein space T 1 ' 1 , the background three- form flux G3 is purely of type 
(2, 1) relative to the complex structure, and the dilaton is constant. We find that the antibrane states 
are dual to solutions where the T 1,1 has been "squashed" (in a way that preserves the full global 
symmetry), the G 3 flux also acquires components of Hodge type (1, 2), and the dilaton is forced to 
run. Additionally, an energy density is generated in the field theory dual consistent with Poincare 
invariance. The squashing, the (1, 2) flux, the dilaton and the energy density are all normalizable 
perturbations of the original theory; they describe a metastable state in the same supersymmetric 
field theory. Using this solution and the holographic dictionary for the KS background determined 
by [IS], we are able to calculate the corresponding expectation values of the relevant operators in 
the SUSY-breaking state of the dual field theory. 

The organization of the paper is as follows. In §2, we warm up by studying linearized supergravity 
equations in a background of the form AdS^ x X 5 . While in a conformal theory one cannot find 
analogous metastable states, many of the features of the solutions to the linearized equations will 
carry over to the confining case. In §3, we present and solve the linearized IIB equations for the 
background of interest, obtaining the coefficients of the linearized perturbations by matching to 
antibrane sources, and describe the features of the solution as well as examining its contribution 
to the force on a probe D3-brane. §4 contains a discussion of the corresponding supersymmetry 
breaking in the dual field theory. We conclude with a discussion of directions for future exploration 
in §5. The appendices contain several calculations whose results we cite in the body of the paper. 

2 Perturbations to Conformal Backgrounds 

Our ultimate interest is in D3-brane-perturbed, supersymmetry-breaking configurations in the non- 
conformal background of the Klebanov-Strassler solution. Before tackling those solutions, however, 
we warm up by considering an analogous computation in the simpler case of an AdS background. 

2.1 Anti-de Sitter backgrounds 

The background geometry around which we perturb is the AdS^ x X 5 solution to type IIB string 
theory, where X 5 is an Sasaki-Einstein space. (In fact any Einstein space will work, although a 
non-Sasaki-Einstein X 5 preserves no supersymmetry.) The nontrivial fields are the metric and RR 
5-form, 

r 2 R 2 ~ Ar 3 

ds 2 = —ri^dx^dx" + —(dr 2 + r 2 dtt 2 x ) , F 5 = (1 + *)— -dr A vol R4 , (1) 
R A r z g s R 

where we defined vo1r4 = dx° A dx l A dx 2 A dx 3 , and R 4 = 47r 4 ^ s A^/vol(A 5 ) with N the quantized 
five-form flux, having set a' = 1. We use the conventions of [17]. The three- form field strengths F 3 
and H 3 vanish, and the axio-dilaton field t = Cq + ie~® is a modulus that can take any constant 
value, with e* = g s the string coupling. For g s small with fixed 't Hooft coupling g s N, string 
corrections are suppressed, and furthermore for g s N itself large we are in the supergravity regime 
with a' corrections negligible. 

The field theory duals to these solutions are conformal field theories (CFTs), the isometry 
group 5*0(4, 2) of the AdS^ space being realized as the four- dimensional conformal group. In known 
examples the flux parameter N is related to the rank of the dual gauge group(s). For example, when 
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X 5 = S 5 the dual is M = 4 super- Yang Mills with SU(N) gauge group, while for X 5 = T 1 ' 1 the dual 
is an J\f = 1 supersymmetric SU(N) x SU (N) theory with bifundamentals and a superpotential. 
We will review this example more extensively when we generalize it to the non-conformal case in 
section [31 

We are interested in perturbing these spacetimes with the introduction of space-filling D3-branes. 
These D3-branes will fill the noncompact directions, and thus four-dimensional Poincare symmetry 
is preserved. A single brane would then sit at a point in both the radial direction r and the angular 
directions on X 5 ; for simplicity, however, we will consider p branes smeared over the compact 
space with a uniform density. This preserves all the isometries of X 5 and hence provides a more 
elementary solution. 

Since the total D3-brane charge associated to the 5-form is related to a parameter in the dual 
field theory, adding D3-branes alone to a given background would result in a change to which field 
theory is the dual, for example from SU(N) to SU(N — p). In order to remain within a fixed 
theory as we perturb, we accompany the p D3-branes with p D3-branes so that the total charge at 
infinity is unchanged. We place these branes at the same location as the antibranes, neglecting in 
the supergravity background the tachyon mediating brane/antibrane annihilation. 

The AdS case has a number of features related to the conformal properties of the dual that 
actually make the solution to this ansatz more subtle to interpret than what we will find for the 
non-conformal case in section [3], despite being technically simpler. The first of these is that we 
cannot obtain static solutions with D3-branes located at any nonzero value of the radial coordinate 
r of the AdS spacetime, because the combination of the background warp factor and 5-form pulls 
them towards r = 0; at best for r > one could expect to find a time-dependent (non-Poincare- 
invariant) solution, with the D3-branes rolling towards r = 0. Thus we cannot include explicit brane 
sources that are static at a finite r = tq, as it would lead to a non-conserved energy- momentum 
tensor and hence to inconsistent Einstein equations. (This is in contrast to the non-conformal case, 
where static D3-branes at a nonzero r will be possible.) 

The alternative, having D3-brane sources at r = 0, is physically admissible but impossible to 
treat directly in perturbation theory in 1/r. Instead, we shall proceed by solving the homogeneous 
fluctuation equations without sources, focusing on the solution far from any branes, and then 
attempting to match to a D3-brane solution by calculating the mass and charge at infinity. A direct 
calculation of the mass will work well for the non-conformal case; for the case at hand we find that 
the mass vanishes, and we will need to look at the limit of branes in flat space to justify the solution. 
This limit will also provide evidence for the identification of the solution as indeed associated to 
antibranes at r = 0, and will generate for us the higher-order corrections in 1/r. 

The other subtlety related to conformal invariance is that in a CFT, the vanishing of the trace 
of the stress tensor: 

T> = 0, (2) 

is an operator equation, and consequently applies to every state in the theory. Any purported 
static, SUSY-breaking antibrane state, on the other hand, preserves Poincare invariance and so any 
VEV for the stress tensor must obey (T^) = trj^ for some t. But then the tracelessness condition 
([2]) requires t = 0; there is no way to preserve Poincare invariance in a CFT while generating an 
energy density. (Of course, dual geometries like a black brane may source the stress tensor; these 
geometries have a horizon and so put time on a different footing, thereby breaking the Poincare 
symmetry and avoiding contradiction with this discussion.) 

This would be a surprise if we expected a state dual to static D3-branes at a nonzero radial 
coordinate r , since one expects such a state to break SUSY at a finite energy scale corresponding to 
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the radial position, and in a globally supersymmetric theory such a breaking must lead to an energy 
density. As we just argued, however, solutions can only correspond to static antibranes if they are 
all the way down the throat at r = 0, corresponding to a vanishing energy scale. We shall indeed 
see that the only solution consistent with our ansatz for the AdS case leads to no energy density 
in the dual, and we will give an interpretation of the solutions as corresponding to brane/antibrane 
pairs at r = 0. 

2.2 Perturbation to the anti-de Sitter background 

Perturbing about this background, we consider an ansatz that preserves the full SO(3, 1) Poincare 
symmetry. Moreover, the branes are smeared over the compact space, thereby preserving the 
isometries of X 5 . For the specific case of X 5 = S 5 this preserves a full 5*0(6) and requires that the 
metric and 5-form ansatz take the form, 



ds 2 = e 2A{r) Vllu dx»dx u + e 2B(r) dr 2 + e 2B(r) r 2 dQ 2 x , F 5 = (1 + *) d^- A dx° A dx 1 A dx 2 A dx 3 , (3) 



where dfl x is the Einstein metric on X 5 . 

For other X 5 the isometry group is smaller and in principle a more general ansatz could be 
considered, but we shall keep this ansatz as the only one that applies to all AdS backgrounds, 
expecting that the smeared antibrane solution will be universal across all choices of X 5 . (We will 
make use of an ansatz that has more freedom to both warp and squash the compact space in 
sectional) 

We also take the 3-form field strengths F 3 and H 3 to vanish, as they would violate the global 
symmetry in the S 5 example, and in any case are not sourced directly by the D3-branes. This in 
turn requires the axion/dilaton, another field not sourced directly by the antibranes, to be harmonic 
according to its field equation, and we choose it to remain constant. (A decoupled fluctuation of 
the dilaton alone also exists, but we will not write it down; a similar mode will show up in the 
non-conformal case.) In principle we could have had different coefficients e 2B , e 2C for the dr 2 and 
dVt 2 x terms, but we used our coordinate freedom to redefine r so as to set B = C, which will simplify 
the calculations. 

We turn now to the field equations. Consider first the 5-form. Its Bianchi identity (equivalent 
to the equation of motion by self-duality) gives 



where the sign of the right- hand- side depends on whether a brane is D3 or D3, y runs over the 
six transverse dimensions to the Minkowski slices, and g is the unwarped metric on the transverse 
six-space ds 2 = dr 2 + r 2 dVL 2 x . We can integrate both sides: 



d r (e- 4A+4B r%a) = £ ±r 5 (2^T 3 ) 



(4) 



3— branes 




(5) 



where N enc is the charge enclosed out to radius r, and so 



n —5 ■ 

o r a — re 



5„4j4-4B 



2n 2 1Q T 3 N e 
vol(X) 



(6) 



Using 2k\ T 3 = g s {2n) A we verify that this is satisfied for an unperturbed AdS solution. 
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Meanwhile, we obtain the following independent Einstein equations. First from the non-compact 
components, 

A" + 4A' 2 + AA'B' + ^A' = ^e~ 8A (d r a) 2 + V 4B (4,T 3 ) £ , (7) 

^ ~ 3— branes v 96 

while from the rr components, 

-4A // -5S /, -4A /2 + 4A / S / -^' = -ie- 8A (^) 2 + ^e- 4B (^ T^ £ , (8) 

^ — 3— branes v 96 



and finally from the compact components, 

B" + AB' 2 + AA'B' + -A' + -B' = ~]e- 8A (d r a) 2 --e- 4B (n 2 w T 3 ) £ (9) 



^ ^ ' 3— branes V yo 

In the Einstein equations both D3-branes and D3-branes contribute with the same sign; we have kept 
the brane sources for completeness, but will ultimately drop them. It is also useful to take a linear 
combination of these equations to remove the A" and B" terms; the combination 4([7]) + + 50 
accomplishes this, and happens to remove the 5-function sources as well. The result is the relation 

12A' 2 + 205' 2 + AOA'B' + —A' + —B' = --e~ SA {d r a) 2 . (10) 

r r 2 

Consider the linearization of the above equations to first order in fluctuations about the background 
AdS solution © e iA ° = e~ AB ° = a = r 4 /R\ We write 

4 

A = \og^ + 6A, B = -log- + 5B, a = ^ + 6a. (11) 

Examine first the integrated five- form equation ([6]). The zeroth order part is satisfied. At first 
order, we obtain 

3 

d r Sa = 16^-(5A-5B). (12) 
R 

We shall use this to eliminate the 5-form contributions from the Einstein equations. From now 
on we drop the delta-function sources, anticipating a solution far from the branes. Linearizing the 
Einstein equations, we find: 

M» + 9^ + 4^ + 32^=0, (13) 

ry ry ry A 

45 A" + 12— + 55B" + — + 32^ = , (14) 

ry ry ry £ 

5B" + 5 5 —- 32 S 4 = 0, (15) 

ry ryZi 

Also useful is the linearized version of the relation fTTUl). which is a linear combination of these three, 



rSA' + -rSB' - -SB = . (16) 
3 3 

Of the three linearized Einstein equations ([TBI . (|14p and (|15|) only two are independent thanks to 
the coordinate invariance of the theory. We choose as our two independent equations the Einstein 
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equation with indices on the compact space (IT51) as well as the relation ffTET) ; it is easy to check they 
imply the other two equations f[T5j) and (fT41) . 

Equation (TT5T) is a linear second-order equation for SB only, and has the obvious power-law 
solutions 

SB = B 4 r 4 + B_ s r' s . (17) 

The relation ( |T6|) then determines the mode <5/L up to a constant: 

<L4 = A - 5 4 r 4 - 2B_ 8 r- 8 , (18) 

and any Aq may be absorbed into a redefinition of the 4D coordinates. Thus f|T7|) . ffT8l) with the 
5-form determined by ffT2l constitute the unique solution to the fluctuation equations around an 
AdS background consistent with the symmetries we have imposed. 

The scalings of r 4 and r~ 8 for this mode are consistent with a fall-off we expect for the non- 
normalizable and normalizable solutions of a field associated to a A = 8 operator Og in the field 
theory dual. Such a mode is known to exist in the KK spectrum of AdS$ x S 5 as well as AdS^ x T ' ; 
indeed it is a universal fluctuation present for all AdS$ x X 5 geometries [18J, owing to the fact that 
it is constant over the compact space and hence is indifferent to the details of the spectrum of 
harmonics. In the M = 4 case, this operator is known to have the form 

O s = Tr {W 2 a Wl) + . . . ~ Tr (F^) 4 + ... (19) 

Thus the leading fluctuation in the metric consistent with the symmetries corresponds to this 
universal A = 8 operator Os- The r 4 mode is non- normalizable, blowing up as the boundary is 
approached, and corresponds to a perturbation of the field theory dual by the Og operator. Hence 
the leading normalizable perturbation to the background AdS$ x S 5 geometry is the r~ 8 mode: 

SB = B r~ 8 , SA = -2B r- 8 . (20) 

This corresponds to an expectation value for the operator G 8 . 

Let us note an operator that is not sourced or given a VEV by this solution: the stress tensor 
T^. As already described, there is no nonzero value of T^ v consistent with both Poincare symmetry 
and conformal invariance, and hence finding zero is what we would expect. Thus there is a vanishing 
energy density in the field theory dual, and hence unbroken supersymmetry, for any value of the 
coefficient of the linearized solution. 

The lack of energy density can also be seen by considering how this perturbation affects the 
appropriate generalization of the ADM mass for anti-de Sitter space. Using the formula from 
Appendix A, we discover that this mass vanishes: 

Eadm oc lim S r~ 4 -> 0, (21) 

r— >oo 

and so this perturbation does not change the total mass of the geometry as compared with the 
background AdS space. In general only a A = 4 operator has the correct scaling to contribute a 
finite amount to the ADM mass. 

Since supersymmetry is not broken and the ADM mass is not shifted, connecting solutions with 
(Os) 7^ to brane/antibrane perturbations may seem like an obscure thing to do. While such an 
interpretation will be much clearer in the case of confining throats discussed next, it is however 
still possible to formally derive our linearized solution from brane/antibrane probe perturbations, 
as follows. 
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Solutions for Poincare-invariant, non-BPS p-branes asymptoting to flat space instead of anti-de 
Sitter space were obtained by Brax, Mandal and Oz [T3]. In Appendix B, we demonstrate that 
the CVperturbed solution found above can also be obtained as the near-horizon geometry of a 
configuration corresponding to N + p D3-branes and p D3-branes all sitting at the origin, in a 
particular scaling limit; the coefficient B is found to go as B ~ p/N, as expected. The overall 
ADM mass, nonzero for the flat space solutions, is scaled to zero in this limit. This derivation also 
produces corrections to the perturbation found here subleading in p/N and at higher order in 1/r. 

Thus it appears that the Os solution is indeed related to brane/antibrane configurations, but 
ones with the perturbing branes sitting at r = 0. We may interpret the motion of the D3-branes all 
the way down the throat as the gravity dual of brane/antibrane annihilation, since the energy scale 
in the dual has retreated to zero, and hence the supersymmetry breaking has disappeared. This 
issue does not occur in the non-conformal case, to which we now turn. 

3 Antibranes in the cascading geometry 

We now consider our primary interest — studying supersymmetry-breaking antibrane perturbations 
of the cascading Klebanov-Strassler geometry associated to fractional branes on the conifold. This 
background differs in two ways from the anti-de Sitter backgrounds. First, conformal invariance is 
broken at all scales (the gravity solution differs from AdS at all values of r) , and consequently there 
is nothing forbidding a Poincare-invariant perturbation that generates a nonzero energy density. 
And secondly, the dual theory is confining with a mass gap, realized in the gravity configuration 
by a smooth termination of the throat at a certain minimal warp factor at r = tq. As a result, a 
D3-brane may achieve a metastable static configuration by relaxing to r , corresponding to SUSY 
breaking at a nonzero but exponentially small scale in the field theory dual. We thus expect 
that asymptotic gravity solutions which are Poincare invariant and carry non-zero energy density 
(determined by the scale of SUSY breaking) will be achievable, and indeed that is what we find. 

In this section we obtain the linearized perturbation to the background associated with the 
addition of brane/antibrane pairs (arguing in §3.2 that this allows us to read off also the asymptotic 
solutions relevant to the antibrane-only states of [7]), and identify the nonzero energy density 
resulting from the supersymmetry breaking. We discuss the geometry of the solution and the 
associated three-form flux, finding that the SUSY breaking results in a "squashing" of the compact 
Einstein space T 1,1 , a running of the dilaton and the generation of a (1, 2) component of the three- 
form flux. We also discuss the fate of this perturbation in the conformal limit, finding consistency 
with the results of §2. 

3.1 The Klebanov-Tseytlin geometry and dual field theory 

The full geometry of the cascade solution is that of Klebanov and Strassler (KS) [6], which consists 
of a warped deformed conifold; the deformation brings the geometry to a smooth end in the infrared, 
which is associated with confinement and a mass gap in the dual field theory. 

We shall focus on solutions at the UV end of the geometry, far from the brane sources which we 
place in their metastable configuration at the bottom of the throat. The UV end of the KS throat 
reduces to the solution of Klebanov and Tseytlin (KT) [15]. Although the KT geometry does not 
include the smooth ending at the tip of the throat, instead continuing down to a singularity at small 
r, it remains non-conformal and must be able to support solutions corresponding to perturbations 
of the true KS throat. Since the full KS geometry can smoothly bring D3-branes to a halt at a 
finite radial distance corresponding to a nonzero energy scale in the dual, it is natural to expect 
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that solutions in the UV exist corresponding to the large r limit of the supergravity fields produced 
by the antibrane configurations; this is indeed with what we find. We will identify the perturbation 
by calculating the associated generalization of an ADM mass. 

The KT background is a deformation of the AdS§ x T 1 ' 1 geometry associated to iV D3-branes 
on the conifold by three-form fluxes corresponding to the addition of M fractional D5-branes. The 
metric is 

ds 2 KT = e^^dxW + e~ 2A ^\dr 2 + r 2 e 2 + r 2 + e%) ) , (22) 



i=i 



where the one-forms associated to the compact space T ' can be taken to be , 

1 2 11 

H = o (<# + /J cos Oid^i) , a dl = —f=dOi , = —?= sin 6^0; , (23) 
3 i=1 V6 V6 

for i — 1,2, and the warp factor is 

_ AA 27ng s ( AT 3g s M 2 l\ 

6 = { + -^(M^Aw) + J , (24) 

where ruv is the scale at which the theory is defined. There is a naked singularity where e~ AA = 0, 
which is simply a relic of the fact that the solution is only valid for sufficiently large distances from 
the tip of the cone; the complete KS solution removes this singularity via deformation. The self-dual 
5-form flux is 

F 5 = (1 + *)F S , F 5 = -277rJV e// (r)vol T i,i , N eff (r) = (n + \og(r/r uv )^j , (25) 

where the volume form on T 1 ' 1 is voI^m = A eg 1 A A ee 2 A e^ 2 . Finally there is 3- form flux, 

9M . , 9g s M . . . , , 

F 3 = — H A {e 01 A - e d2 A e^J , £ 2 = (e 0l A - e 02 A e fe ) log(r/rav) ■ (26) 

The axio-dilaton r = Cq + ie~* remains constant. 

The preserved symmetries of this background are 5*0(3, 1) of the four- dimensional along with 
the S77(2)i x SU{2) 2 x Z 2M x Z 2 of the T 1 ' 1 ; the SU(2)s act on the two-spheres associated to 6*i, 4>i 
and 2 ; 02, and the Z 2 exchanges them. The metric possesses a U(l)^ acting on the -^-coordinate, 
which is respected by the other field strengths but broken to Z 2 m by an Y concrete expression for 
the potential C 2 [20]. This is broken further to Z 2 in the full KS solution, reflecting the fermion 
condensate in the IR. 

In the M — > limit we revert to AdS§ x T 1,1 , corresponding to N D3-branes sitting at the tip of 
the conifold. The dual field theory to this conformal limit is an M = 1 superconformal field theory, 
which can be realized as an SU(N) x SU (N) M = 1 gauge theory with bifundamental matter 
and Bj in the (N, N) and (N, N) of the gauge group. SU{2)\ x SU{2) 2 is a global symmetry with 
Ai and Bj transforming as doublets of the first and second factor respectively, and moreover there 
is a superpotential 

W = Xe lk e jl Tr A i B j A k B l . (27) 

U(l)^ is restored in the conformal limit as the R-symmetry. 

This superconformal theory has two complex moduli parameterizing the space of marginal de- 
formations, which may be thought of as gauge couplings and theta angles for each SU(N) factor; for 
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given values of these couplings, the superpotential parameter A is adjusted to reach the fixed point, 
where Ai and Bj have dimension 3/4 and the superpotential is marginal. The complex moduli are 
dual on the gravity side to the axio-dilaton and to /^(.Ba +2C2), respectively, and their associated 
operators, along with the energy-momentum tensor dual to the four- dimensional metric, make 
up the exactly marginal (A = 4) operators of the theory. 

For nonzero M, the theory becomes an SU(N) x SU (N + M) gauge theory, with Ai and Bj still 
bifundamentals of the new gauge groups and with the same superpotential. Conformal symmetry 
is broken, the R-symmetry U(l)^ is broken to Z2M by anomalies, and one of the two moduli is 
lifted. The simple way to think of this is that the sum of the gauge couplings l/gf + l/g^ remains a 
modulus dual to the constant value of the dilaton r, while the difference begins to flow with energy 
scale, 

-a - \ ~ — / B 2 ~M \og(r/r uv ) . (28) 
9i 92 9s Js* 

We shall refer to the operator associated to the sum of the gauge couplings as + , and the one 
associated with the difference as CL. As we discuss later, properly both couplings also involve the 
superpotential parameter A. 

The RG evolution exhibits the beautiful phenomenon of self-similar flow called the duality 
cascade, where the effective number of colors N e ff diminishes as the scale decreases, corresponding 
to successive Seiberg dualities that progressively reduce the ranks of the gauge group SU(N) x 
SU(N + M) — > SU(N — M) x SU(N) and so on, while preserving the structure of the theory at 
each step. 

3.2 Relation to metastable antibrane states 

Before proceeding, we should note the relation of these states to the states of [7j. There, one 
considered a warped throat to be part of a compact geometry, with 

[ F 3 = M, { H 3 = -K. (29) 

J A JB 

Adding p <C M D3-branes to this geometry was then argued to result in metastable supersymmetry- 
breaking vacua. How does this compare to the case we study with p ordinary D3-branes as well as 
the p D3-branes? 

Consider the issue of total charges. Because of the Bianchi identity for F 5 , these fluxes contribute 
MK units to the total D3-brane charge Qt>3, which schematically reads: 

Q D3 = J dh = N m - Afe + J H 3 AF 3 = N m -N m + MK, (30) 

The tunneling decay mechanism of the p D3-branes must conserve this total charge. This proceeds 
by the dynamical nucleation of an NS5-brane bubble which wraps the A-cycle; this brane acts as 
a source for H 3 and shifts K — > K — 1. This alone changes the total charge (!30|) by —M, but 
the process compensates by also creating M D3-branes ending on the bubble (and filling the space 
"outside"), p of these newly-created D3-branes annihilate the antibranes, leaving a supersymmetric 
configuration with M — p D3-branes, background flux giving rise to charge — M, and no D3- 
branes. These together have D3 charge N — p, while the F 3 flux remains M. Consequently, the 
p-D3-brane state should be viewed as a metastable state, not in an SU(N + M) x SU (N) theory, 
but in the SU(N + M — p) x SU (N — p) supersymmetric gauge theory associated to N — p D3s 
and M fractional D5s at the conifold. 
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Imagine now a noncompact case where we let N + p be an integer multiple of M, instead of 
N = KM. The asymptotic geometry we use to study this theory, after the addition of p probe 
D3s and p probe D3s to the background characterized by iV (but forbidding the brane/antibrane 
tachyons to condense), should have identical conserved charges compared to the theory of N + p 
D3s and M fractional 5-branes with p D3-brane probes. Since now iV + p is divisible by M, we 
are back in the situation considered in [TJ, and the asymptotic solution we find should agree with 
that for p probe D3-branes in the confining Klebanov-Strassler geometry as far as the asymptotic 
charges are concerned. Therefore, if we find a solution with the correct values of the asymptotic 
charges, motivated by adding p D3/D3 pairs to the geometry with N D3s and M fractional D5s, 
we can be confident it coincides with the supergravity solution for the metastable D3-brane states. 

The reader may have found the preceding argument somewhat confusing; let us give a slightly 
different explanation of the same fact. We could formally find a supersymmetric solution with zero 
vacuum energy, by adding p "imaginary branes" with negative D3 charge and negative D3 tension 
to the geometry with N + p branes; it would be supersymmetric because these "imaginary branes" 
satisfy a BPS condition relating their charge and tension equivalent to the relation satisfied by the 
N + p "real" branes. The solution we are interested in differs from this one, by adding p D3-branes 
and p D3-branes. (The p D3s can then cancel against the p "imaginary branes" as sources for charge 
and energy density). Hence, our solution of interest should have an energy density characteristic of 
2p (warped) D3s, and the same charge as the solution with N units of D3 charge. This argument 
was first given by Maldacena and Nastase in a rather different context (supersymmetry breaking in 
the gravity dual of a 3d M = 1 supersymmetric Chern-Simons theory) [21 j . 

3.3 Perturbation ansatz in KT background 

We now wish to consider the response of the cascading geometry when it is perturbed by a D3- 
brane/D3-brane source. 

One way to approach this problem would be to explicitly place branes and antibranes at the 
base of the full KS geometry and discover the response of the system in the UV. This would entail, 
however, somewhat tedious matching between the IR and UV regions. Instead we shall pursue the 
same philosophy as in the AdS case, and consider sourceless fluctuations in the UV regions — that 
is, in the KT geometry — and then match these to an ADM mass to verify that they correspond to 
the addition of tension without charge. We shall also demonstrate explicitly that branes will source 
the particular fields that are turned on. 

As in the conformal case of the previous section, we will simplify matters by smearing the 
perturbing branes over the compact space such that the full symmetry of T 1,1 is maintained. We 
thus preserve Poincare invariance as well as SU{2)\ x SU(2) 2 x Z 2 , and we keep the metric U(l)^ 
symmetric as well. Our ansatz for the metric is thus constrained to be 

ds 2 = r 2 e 2a(r) v&fl(f/ + e -2a(r) ^ + ^(^2 + ^(r) jj^ + q\ ? ^ 

where we have found it convenient to extract a factor of r 2 from the warp factor. Additionally, we 
have used our freedom to perform coordinate redefinitions of r to relate the coefficient of the dr 2 
term with the coefficient of the Minkowski part. 

Notice that even maintaining all symmetries due to the smeared branes, it is possible to introduce 
a new metric function, the ratio e b ~ c between the length scale of the circle acted on by U(l)$ and 
that of the two S 2 s. We shall find this degree of freedom to be essential in what follows. 
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For notational convenience we define the following rescaled constants: 



N = 27nN, M=-M, (32) 
2 

and for the remaining fields we take, 

F 3 = Me^A(e Sl Ae^-e fe Ae^), B 2 — k(r)(eg 1 A — eg 2 Ae^), (33) 

F 5 = dC A -C 2 AH 3) C 4 = ^UoW, (34) 

9s 

$ = $( r ), C = 0. (35) 

where again vol^i = dx° A da; 1 A dx 2 A dx 3 . For the five-form, we use the forms of the potentials 
and the Bianchi identity to obtain, 

F 5 = (1 + *) T h , F 5 = -C 2 AH 3 = -N eff (r) vo1 t m , N eff (r) =N + 2Mk(r) , (36) 

and the potential C 4 is determined by self-duality to be 

d r a = -g s N eff (r)r 3 e 8a - b - 4c . (37) 

We now look for solutions to this ansatz. The field equations of Type IIB SUGRA are [T7] , 

V 2 $ = e ^d M C d M C - 9 ^(H MNP H MNP -e^F MNP F MNP ) (38) 

V M (e 2 *d M C ) = -^H MNP F MNP , d*F 5 = -F 3 AH 3 , (39) 

o 

d*(e*F 3 ) = g s F 5 AH 3 , d * (e~' s 'H 3 — C e®F 3 ) = —g s F 5 A F 3 , (40) 

I e 2 * g 2 ~ ~ p 
Rmn = -9m^9n^ H — —OmCqS Co + —FmpqrsF n 

II yo 

+ ^(^/PQ^/ Q + e 2$ F MPQ F/ Q ) (41) 
- g ^9MN(H PQR H^ R + e^F PQR F p ^), 

where M, N = 1, 10 and F 3 = F 3 — C$H 3 . Checking the consistency of our ansatz, we see that 
the axion Cq can indeed be set to zero as long as F 3 and H 3 don't have any components which are 
polarized in the same directions. We have already satisfied the F^ Bianchi identity, and the F 3 and 
H 3 Bianchi identities are trivially satisfied by our ansatz as well. The F 3 equation of motion can 
be seen to be be satisfied using F 5 A H 3 = and 

* e"F 3 = Mr 3 e* +4a " fe vol R 4 A dr A (e^ A e 6l - e 02 A e 02 ) , (42) 

which is closed. 

The equations that remain are the dilaton, H 3 , and Einstein equations. The first two take the 
explicit form 

^* = "7^(^ " 7^) ' (43) 
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l^^ M )=S,M^(J^), (44) 

In addition, we have three distinct Einstein equations, which are 

p 9y Nl ff g s e-* 4 2 r 2 e 2 *M 2 2 P n\ r 2 T 3 5{r - r ) 

xx 4 g-12a+26+8c 4g-8a+4c V g26 ' " t " 2e~ 2a y/Cfe ' 

for the component, 

p _^ 2 ^// 3e 2 *M 2 2p^ o r 3 e 2 ^(r - r ) 

^ 36e- 8 «+ 8c 36e" 4 «+ 4cl e 26 ) 18e 2 ^ ' 1 J 

for the ipip component, and 

p gj ^// ^ e 2 ^M 2 2^ 2 T 3 e 2 ^(r - r ) 

ee 24e- 8a+2 ''+ 6c 24e" 4a+2cl e 26 J V2e 2a Jg^ ' 1 J 



for either 9i9i component; the full expression for the components of the Ricci tensor for the metric 
( !3T|) can be found in Appendix C. 

We can gain some intuition concerning the way the various equations couple with one another. 
One may see that the RHS of the xx-Einstein equation (145]) acts as a source for e 4a . Substituting 
this result into the other two Einstein equations ()4*6l) . (I4T1) . one can show that e 2c is not sourced 
and it is consistent to set c = 0. 

However, in general the function b is nontrivial. This is interesting because e b ~ c represents the 
relative warping between two factors of T 1 ' 1 ; the four dimensional space comprised by the S' 2 s on 
the one hand, and the ^-circle acted on by U(l)^ on the other. Thus if c is trivial, but b is not, 
then the compact space is squashed. 

In fact, we find that the squashing e 2b is sourced by the combination of fluxes, 

^ k ' 2 ~ -^-)- (48) 

Observe that the source for the squashing term ( )48|) contains the same function as the right-hand- 
side of the dilaton equation ( )43|) . which referring to ( 1381) can be shown to be proportional to the 
flux combination 

F 3 ■ F 3 - e- 2 *H 3 ■ H 3 oc Re G + ■ G_ , (49) 

where as usual we have defined the complex flux G 3 and its imaginary-self-dual and imaginary- anti- 
self-dual combinations: 

G 3 = F 3 - tH 3 , G± = iG 3 ± H G 3 . (50) 

The KT background has flux that is purely imaginary-self-dual, G_ = 0. Thus we find that 
a perturbation that generates an imaginary anti-self-dual flux in addition must both source the 
dilaton, and squash the compact space. 

We shall find that exactly this occurs. The antibrane itself does not directly source the flux, but 
sources the metric warp factor e 2a . This in turn acts as a source in fj4*4l) and produces a nontrivial 
perturbation to k corresponding to flux that is no longer imaginary-self-dual. Finally, this flux 
causes the dilaton to run and the T 1,1 to squash. We thus find the interesting result that even 
branes smeared so as to break no symmetries of T 1,1 in the UV still end up modifying its shape. 
The compact space is highly symmetric, but not symmetric enough to prevent this. 
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3.4 Linearized solution and its properties 



We now explicitly solve equations (l4~3l - (l4"Tj) to linear order. To do so, we expand the functions, 
e 2a , e 2b , e 2c , $, k, in a power series around infinity at 1/r = 0. We consider an ansatz where each 
function has corrections at order 1/r 4 as well as log r/r 4 ; higher order terms are neglected and the 
coefficient of the perturbation is kept only to linear order. We find a consistent solution to the 
linearized equations depending on two parameters, which we call S and 0: 



-la 



,26 



-,2c 



\g s N + ~(g s M) 2 + -(<? s M) 2 logr 
1 



32 

! + ^, 
1, 

g s M\ogr + 

l0g# S + ~7 [< 



1 - 13 1 

-g s N + — (g s M) 2 + -ig s M) 2 logr ) S - i~{g s M) 



64 



16 




^g s M+^g s Mlogr\s 



;9sM< 



3 S log r] 



(51) 

(52) 
(53) 
(54) 

(55) 



We claim that this is in fact the unique solution consistent with our ansatz. We provide additional 
justification for this by identifying this solution as a particular member of a class of solutions 
generated by Aharony, Buchel and Yarom [16] . which are claimed to exhaust the solutions to 
a class of perturbations around KT broad enough to include our ansatz; the comparison to the 
notation of [16] is made in Appendix D. 

We would like to calculate the total mass of this solution. Since the cascade geometry is asymp- 
totic neither to flat space nor to AdS, we cannot use ordinary expressions for the mass specialized 
to those geometries. In Appendix A we review the Hawking-Horowitz mass [22J which can be com- 
puted in more general spaces. We find for the contribution of the perturbed background to the 
normalized extrinsic curvature integral, 



Af / K 



3r 4 + 



3 4a 



D 8ao 



G? s m)V 



D 4a 



-s 



—g s N+—(g s M) 2 + -{g s Mf logr )S - ~{g s M) 



16 

{g s Mf 



32 



(56) 



13 



-g s N + -(g s M) 2 + -(g s M) 2 logr S - -(g s M) 



32 



64 



4 



16 



where we have suppressed a volume factor V^vo^T 1 ' 1 ) 
factor, 



and cio is the background value of the warp 



-4a 



\g s N + \{g s Mf + \{g s M) 2 logr . 



(57) 



Subtracting the background value M J /Co for the unperturbed cascade removes the first two terms 
which are divergent at large r. We then take the r — > oo limit; the warp factor e _4a ° diverges as 
logr. Most terms vanish, and we find the finite mass (restoring volume factors) 



E 



1 



ItiG 
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M j{K - /C ) = V 3 voKT 1 - 1 ) S 

J 167rGio 



(5c 
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where Giq = k1 /8tt is Newton's constant in ten dimensions. We see that the mass is directly 
proportional to the S perturbation, while the perturbation is not involved at all. Meanwhile, the 
value of the D3 charge is trivial due to the definition of F§; it is simply N e ff(r). All corrections to 
k(r) vanish as a power 1/r 4 in the r — > oo limit, so this solution does not change the asymptotic 
charge of the cascade. 

We identify the coefficient S as the correct parameter associated to the brane/antibrane per- 
turbation. The mode can be seen to reduce in the AdS limit (M — > 0) to the usual independent 
fluctuation of the dilaton, related to the ability to turn on a vacuum expectation value for the dual 
operator + in the absence of 3-form fluxes. As the brane/antibrane probes do not source the 
dilaton directly, we set this parameter to vanish, as we implicitly did in the AdS case: = 0. Note, 
however, that although this means there is no non-logarithmic perturbation to the dilaton field, 
we are still left with a dilaton perturbation proportional to (log r/r 4 ) S. This is because, as was 
emphasized previously, the perturbation sources IASD flux, which in turn squashes the T 1 ' 1 and 
forces the dilaton to run. This logarithmic perturbation is that running. 

The numerical value of S is determined by the boundary conditions on the fields at the IR tip 
of the geometry. Although we do not solve this precisely because we choose to avoid the details of 
the IR part of the geometry, we may still use the way a source contributes to the field equations 
to determine the dependence of S on the various parameters. The Einstein equation (j45p in the 
presence of the p D3/D3 sources localized at r = r gives 

d r d r e~ 4a + ... = 4^ T 3 r 4 (5(r ~ ro) , (59) 

where = dr 2 +r 2 ds^ ltl is the unwarped metric on the transverse six-space and we have only written 
the second derivative term on the warp factor. Multiplying both sides by 1/r 5 and integrating over 
the transverse six space, we find S ~ (p/N) r 4 @ Thus, S is warped by r 4 and proportional to the 
number p of brane-antibrane pairs. 

This scaling of the perturbation S enables us to understand how it behaves in the AdS limit of 
M — > 0. Examining the form of the solution for k(r) (1541 . it naively seems to diverge in this limit 
due to the term proportional to N/M. However, this divergence assumes that S is kept fixed in 
that limit; in fact it is not. The IR scale of the throat r depends exponentially on the fluxes as 
r ~ exp(— 2nN/3g s M 2 ), thus receding to zero in the AdS limit, and likewise 

sending the total perturbation to zero as M — > 0. This confirms explicitly that the supersymmetry 
breaking is exponentially small relative to the fundamental scale, corresponding instead to the IR 
scale of the bottom of the throat. 

This behavior reflects the fact that although the KT solution does not directly include the smooth 
IR tip, nonetheless the KS deformation is the only way to resolve it in a non-singular fashion, and 
hence perturbations around KT "know" that the IR tip is their only consistent resolution. The 
geometry is aware that it cannot support a solution of the kind we have found here once the AdS 
limit is taken — indeed we saw no trace of it in the analysis of section [2] — and the fact that it 
must be scaled to zero in that limit corresponds to the retreat of the confining scale down to zero. 

2 Strictly speaking the quantity in the denominator of p is some combination of the two background contributions 
to the running number of colors, N and g s M 2 . 
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3.5 Brane/antibrane force 

Another check on our linearized solution is provided by a comparison to the calculation in [23] of 
the force between a well-separated brane/antibrane pair in a warped throat. (In |23j the calculation 
was approximated assuming a background AdS geometry; the true conifold case was considered in 
[24|). Placing the antibrane at the cutoff tip at r and the brane at r, for r ^> r , the attractive 
potential between the two was found to be 

V(r) = - 2Tsr8 ° . (61) 
1 J NR*r 4 K ' 

Also, observe the warping proportional to r®. This potential results in an attractive force F{r) of 
magnitude 

To compare with this result, we place a probe D3-brane in the background of our linearized solution, 
which already includes the effects of p D3-branes at r = r , and investigate the resulting force. To 
lowest order in a', the brane probe is governed by a Lagrangian consisting of the DBI and CS terms 



Sd3 — ~T-.i 



J d 4 x r 4 e 4a W 1 + ( -^d il rd»r + T 3 Jd 4 xa, (63) 



where \i = 0, 3. 

Expanding the action, we find in fact that there is zero force on the probe brane at the order to 
which we have worked. In fact this is the result of a precise cancellation; if any of the coefficients in 
the perturbation (!5T!) - (!55l) are tuned away from their actual precise numerical values (which involve 
numbers like 13/64 and 11/16), then a force would result of order F(r) ~ r^/R^r, which falls 
off much more slowly than (j62j) and would have dominated over that result at large r. Only this 
cancellation prevents the cascade geometry from having a much stronger flux- aided brane/antibrane 
force than the result of [23J. 

Given then that the expected force on a probe brane is not the result of the leading dimension 
four operator perturbations that we found, one may ask whether this force exists and what produces 
it. Recall that we have only worked to order 1/r 4 , and hence perturbations vanishing faster at large 
r will be invisible. If one postulates the existence of additional terms at order 1/r™ with n > 4 and 
coefficient T (including possible log terms), but keeps both S and T only to linear order, then one 
must have T ~ p and a force is generated, 

\F(r) ~ T 3 5— — , (64) 

r n ~ 3 N^ ff (r) v ' 

for some constants 5, e fixed by the equations of motion. The correct force law is seen to arise for 
n = 8. 

There are two modes consistent with the symmetries of our ansatz corresponding to higher- 
dimension operators that could be turned on in the background at higher order: a dimension six 
operator C 6 associated with the warp factors e b and e c , and the familiar O s operator discussed 
in section [2j Given the matching of the force calculation with n = 8, it is natural to surmise 
that the operator Og is turned on in the background as well at order p, and is responsible for the 
brane/antibrane force; it should be visible if one were to include perturbations out to 1/r 8 . The 
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operator Oq would then either be absent, or would participate in another nontrivial cancellation to 
avoid producing a 1/r 3 contribution to the force (IMl) . (Note that the S perturbation would also 
correct the 1/r 5 force law at order p 2 .) 

Due to the sensitivity of the 1/r force cancellation to the precise coefficients of the solution, the 
above calculation is a highly nontrivial check on our result; although the dimension eight operator is 
responsible for the force on a probe brane, we see this is only possible with the precise cancellation 
that occurs with terms associated to the dimension four operator. 



3.6 Three-form flux type 



As mentioned before, the KT (and KS) geometries possess imaginary-self dual flux G + , with G_ = 0. 
Moreover, this flux is known to be of type (2, 1) relative to the complex structure of the geometry; 
in general ISD flux can be (2, 1) or (0,3) while IASD flux is (3,0) or (1,2). In this subsection we 
determine the Hodge or index type of the flux perturbation. 
The complex flux G3 takes the form 



G, 



(Me, 



ie 



-*d r kdr)A — 



(65) 



where u 2 = S(eg 1 A — eg 2 A e^). Tensors associated to different Hodge types in this background 
were classified by Herzog, Klebanov and Ouyang |25j. They provided a dictionary to translate the 
forms dr, and 002 we use into the four complex coordinates Zj defining the deformed conifold via 

E 



A. 



dr 

h 



cu 2 



Zi dzi 

3 r d 



dr 
r 



Zi dz{ 

3 r 6 



tijkiZiZjdzk A dz\ 



in terms of which the holomorphic structure is manifest; note the scaling \z\ 
we have, 



(66) 
(67) 

r 3 . Then for (1651) 



° 3 ~ 9^ 



(m + 



T'drp k ) dz, 



(m 



£ijM z iZjdzk A dz\ . 



(68) 



We readily see that the only flux types are a (2, 1) part (hence ISD) proportional to [M + e~®rd r k), 
and a (1,2) part (hence IASD) proportional to (M — e~®rd r k). The product of these quantities 
is the product of ISD and IASD fluxes, and indeed is proportional to the flux combination f|48l) 
sourcing the dilaton and the squashing, which we previously identified as Re G + G_. 

Since the M piece of each G + and G_ is part of the background, the contribution of the 
fluctuation comes entirely from the e~^rd r k terms, and hence the brane backreaction contributes 
fluctuations of (2, 1) and (1, 2) fluxes with equal magnitude. For completeness, the fluxes are 



G+,o 
G-,o 



2i 
979 

2i S 
9r^ r 1 



S f3 N 



2M - — - 



3 N 



2 ai M + i a + Ufl ° Sr , 



(j ]f + -^M + 3Mlogr 1 (e m ZiZjZ m dz k A dz t A dz m ) . 



(etjkiZiZjZradzk A dz x A dz m ) , (69) 

(70) 



We emphasize that the above represents the decomposition of the perturbed flux into the Hodge 
types of the background geometry, denoted by the subscript. However, the squashing factor e b will 



16 



also modify the metric, and so a natural question is whether the perturbed flux is ISD with respect 
to the perturbed geometry. One may show that the fluxes take the form, 



M ± e b e-*rd r k 



ie^±e —I Ay, 



(71) 



in the modified background, which evaluates to 



G + 



S ( 3 N 3- - 

2M-— ^ + -M + 3Mlogr 

r 4 \2g s M 4 6 , 

S /3 N 3 - - \ / 
-j ( --^ + -M + 3M logr 



^ + e- 7 Ay 



r 



2g s M 




(72) 
(73) 



and hence we see the perturbed flux is not ISD (or IASD) with respect to the perturbed geometry, 
either. 

We note in passing that the leading corrections to the KT solution coming from the full KS 
solution, indicative of chiral symmetry breaking, actually fall off less quickly than our perturbation 
far from the tip, as G3 ~ M(ro/r) 3 [26]. Hence our perturbation is not the leading correction to the 
KT geometry; it is, however, the leading term proportional to p and the leading term that breaks 
supersymmetry. An antibrane solution valid over the entire KS geometry should contain corrections 
to the chiral symmetry breaking structure proportional to p as well. 



4 SUSY Breaking in the Field Theory Dual 

Because the metastable D3-brane states have not proved amenable to a direct weak 't Hooft cou- 
pling field theory analysis (for work in this direction in slightly more general geometries, see |llj). 
the characteristics of the supersymmetry breaking that we can learn from the operator vacuum 
expectation values using the AdS/CFT dictionary may be a useful source of insight. 

Even before looking at the VEVs, we can immediately see one interesting feature. As explained 
in [20], the anomaly-related R-symmetry breaking in the Klebanov-Strassler solution from U(l) to 
Z 2 m is visible through the behavior of the C2 Ramond-Ramond gauge field. In our solution, F 3 and 
hence C2 are completely unchanged from their background value. So the SUSY breaking does not 
(at least at leading order) give rise to R-breaking independent from that already present in the KT 
and KS solutions. 

Because the SUSY-breaking does not provide another independent order parameter for R- 
breaking, we conclude that it is either D-term breaking, or F-term breaking which preserves R. 
The reader should note that while 0(1) R-breaking at the SUSY-breaking scale would imply F- 
term breaking (or at least a significant F-component to the breaking), the absence of R-breaking 
is consistent with either D or F-term breaking — with the ISS models [27] providing a canonical 
example where F-term breaking preserves R. 

To learn more about the supersymmetry breaking, let us try to characterize the state in the 
field theory more precisely. Being a set of normalizable fluctuations of the supergravity fields, our 
solution implies a number of vacuum expectation values for the dual operators. The techniques 
of holographic renormalization [28] for calculating expectation values in an asymptotically anti-de 
Sitter background have been generalized to the cascade geometry by Aharony, Buchel and Yarom 
[T6] and our solution can be treated by their techniques; we match our notation to theirs in Appendix 
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D. Using their formulae, we find the expectation values (including the (f) mode for completeness): 



(o + ) 

(0_> 



-2,S - 40 



(74) 
(75) 
(76) 




where 0+ represents the operator dual to e while O- represents the operator associated to 
As with the perturbed solution itself, the vev seems singular in the M — > limit at 
first glance, but actually goes smoothly to zero since S vanishes exponentially ([60]) . 

If one defines the 5D boundary stress tensor as Tf£ = T^/SttG^ with G 5 = Gio/voliT 1 ' 1 ) as in 
[To] , one sees immediately that (Tq®) exactly matches the calculation of the mass density of the 
spacetime presented in (1581) . The nonzero value of the energy density is confirmation that the state 
we are studying indeed breaks supersymmetry, and the scaling ([60]) of S confirms that it is broken 
at an exponentially small scale. 

The running of the coupling associated to O- (128]) generates a quantum violation of conformal 
invariance proportional to the beta function: 



where contributions proportional to the conformal anomaly vanish for a Minkowski background. 
This operator relation must hold in any state, in particular our supersymmetry-breaking state. 
It was verified by [H] that for the beta function at hand, which becomes (3- = M for us0 the 
corresponding relation 



is indeed satisfied, as can be confirmed from (1741) and f!76|) . Note that the dependence of (0 + ) on </>, 
the parameter we set to zero for the D3-brane solution, but which indicates another set of possible 
states in the theory, rules out the appearance of + in such a relation. 

Now, let us return to the question of whether the supersymmetry breaking can be characterized 
as D-term or F-term breaking. For our Poincare-invariant state, the trace of the stress tensor is 
proportional to the energy density; thus (!7H|) is a relation between the energy of the state and the 
expectation value of O-. Hence if we can determine the Lorentz-invariant terms in O- that can 
acquire a VEV, we can characterize the nature of the supersymmetry breaking. Additionally, a 
super symmetric partner of the operator relation (|77l) relates the fermionic partner of O-, call it 
to the supercurrent, which when supersymmetry is spontaneously broken is linear in the goldstino; 
thus we expect to be linear in the goldstino when its auxiliary fields are given their VEV. 

For these reasons it is interesting to understand (9_ better. We have mentioned previously that 
+ and CL can be identified with the sum and difference of the coupling constants of the two 
factors in the gauge group. Properly speaking however, the story is a little more complicated [HE]. 
In the conformal limit, there are three SU(2)i x 5 , f/(2) 2 -preserving dimension four scalar operators: 
the F-terms of Ti^Wy, Tr(W^ 2 ) anc i Tr (A i B k A j B l e ij e kl ), the first two associated to the two 

3 The imaginary parts of 0± operators, dual to modes of Co and J C2, are not given an expectation value in our 
solution. 

4 Note we have denned C_ as directly dual to e - *f?2 with no numerical factors, which leads to the given normal- 
ization. 



(77) 



(TP = -\m(0.) , 



(78) 
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gauge couplings g\ and g 2 and the last one being the superpotential term (j2"7|) . which has associated 
coupling A. There are, however, only two combinations of the three couplings that are invariant 
under all the anomalous global symmetries in the CFT, which we can take as J = (A 2 AiA2) 7V and 
/ = (A 2 /Ai) N where Ai and A2 are the usual complexified dynamically-generated scales for the 
two gauge factors; it is easy to see that J and / do involve the sum and differences of the gauge 
couplings 1/gf and l/g%, respectively. 

The freedom to choose these two parameters in the CFT corresponds to the two moduli of the 
theory; put differently, there is a two-complex-dimensional space of conformal fixed points in the 
three-dimensional space of Ai, A2 and A. These can be identified with the moduli r ee Co -Me - * and 
<p = e~* f S 2(B 2 +iC 2 ) of the gravity dual by matching the action of the Z 2 center of the SL{2, Z) S- 
duality of type IIB string theory, which on the field theory side acts as charge conjugation combined 
with an exchange of the two gauge factors and thus leaves A and B invariant [4J. This symmetry 
preserves both r and J, which we thus match together (calling the associated operator + ), while 
flipping the sign of ip and J, which we thus match (calling the associated operator OJ). More 
strictly speaking, r and ip are associated to log J and log J, respectively. A more detailed discussion 
of the identification of the Z 2 odd modulus with J, using the fact that one can obtain the KW 
theory from a mass-perturbed Af = 2 super conformal field theory, appears in 

Away from the conformal limit, the generalizations of J and / are 



a N—2M 

r 1 V*- 



J EE X ™+M A 3M + N A N-2M ? j = \3M X JV__ V<] ™ . (79) 

For nonzero M, we see I is no longer invariant built solely out of coupling constants, but instead 
becomes field-dependent. This reflects the fact that it is no longer associated to a modulus; the 
difference of the gauge couplings starts to run as in ( 1281) . 

From all this we learn that in the conformal limit, 0_ (unlike + ) does not involve the super- 
potential, and has the simple form 



CL ~ Tr 



(F( + Ai<9Ai + Df) - (F 2 + X 2 d\ 2 + D%) , (80) 



where we have schematically indicated the gauge fields, gauginos, and auxiliary D fields for the two 
gauge factors. Hence a Lorentz-invariant expectation value for this operator must involve solely the 
auxiliary fields in the gauge multiplet, and thus constitutes D-term supersymmetry breaking. 

However, our metastable states only arise in the non-conformal theories with M^O — there are 
no Lorentz-invariant metastable states in the conformal theory. Away from the conformal limit, we 
see from (179|) that the superpotential does begin to mix into CL. Thus in principle, the auxiliary 
fields Fa and Fb may begin to contribute to an F-term supersymmetry breaking component as 
well. While the fact that in the conformal limit (9_ is dominated by D-terms is suggestive of 
D-term supersymmetry breaking, the magnitude of the energy density is small enough that the 
sub dominant F-term mixing could parametrically account for the full result. Therefore, while our 
results suggest that the breaking may be dominated by D-terms, we cannot give a precise argument 
that this must be so. In fact, we note that in a strongly coupled gauge theory where one cannot 
identify the Goldstino uniquely as a component of some (weakly coupled) gauge or chiral multiplet, 
there may well be no invariant distinction between D- and F-term supersymmetry breaking^] 



5 In a theory with a U(l) factor in the gauge group and a Fayet-Iliopoulos D-term, this situation may be different. 
In coupling to supergravity, the FI term causes shifts to gauge charges and may thus distinguish the two possibilities. 
However, our theory is an SU(N) x SU(N + M) gauge theory with no U(l) factors in the gauge group. We thank 
N. Seiberg for discussions of this possibility. 
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5 Conclusion 



Motivated by the desire to find gravity duals of vacua with exponentially small supersymmetry 
breaking scale, we studied linearized perturbations of supergravity backgrounds. For backgrounds 
dual to superconformal field theories, we found a universal dimension eight operator, and a possible 
identification between backgrounds perturbed by a VEV for this operator and brane/antibrane 
configurations, but nonetheless supersymmetry remained unbroken. 

Our primary interest was to study perturbations of backgrounds dual to confining theories, 
concentrating on the KS solution. We identified linearized perturbations corresponding to the 
metastable SUSY breaking states of [TJ. Let us summarize the effects of the SUSY breaking D3- 
branes on the background geometry. The D3-branes directly source the warp factor; in fact, they 
require at least two warp factors (or, an overall warp factor and a "squashing" factor). The warp 
factor then acts as a source for the background three-form flux, and we find that both G± are 
activated. The combination of both G + and GL flux then sources both the dilaton and the squashing 
factor. As a result, the final solution differs in several precise ways from the supersymmetric 
solutions: the T 1,1 is squashed, the three-form flux includes (1,2) components, and the dilaton is 
forced to run. 

There are a number of promising directions for future work. Soft-breaking terms induced by 
SUSY-breaking three- form fluxes on D3-brane probe fields were found in [30J . These works assumed 
a single warp factor, i.e. a conformally Calabi-Yau metric ansatz (and as a consequence, worked 
out soft terms in flux backgrounds that in general do not solve the classical closed string equations 
of motion). Therefore, their results cannot be literally applied to our solutions. It would be 
interesting to modify these analyses of the soft terms to make them applicable to more general 
backgrounds like the ones we found. It has been suggested that antibrane SUSY breaking in 
warped throats is naturally sequestered, with detailed arguments appearing in [31]. If this is so, 
one would expect the leading approximation to the soft terms to vanish in many circumstances. 
Additionally, the background we describe here and related configurations may have applications in 
phenomenologically-motivated settings like that studied in [32] . 

Compactification of supergravity solutions of the form presented here is still a nontrivial open 
problem. The simplest IIB backgrounds in a compact setting have a conformally Calabi-Yau metric 
and no (1,2) fluxes [33]. It may be possible to find more general compact, classical backgrounds 
which allow more general supersymmetry-breaking possibilities such as the one presented here. 

There are also several clear steps one could take to make our solutions more complete: 

• We worked to leading nontrivial order in 1/r. One could work to higher orders in 1/r, to e.g. 
capture the leading force on a brane probe at large r, which was approximated in [23J (recall that 
we found the expected "miraculous" force cancellation at the order we worked). Similarly one could 
also keep terms of higher order in the p/N expansion. 

• In order to simplify the analysis, we "smeared" the D3-branes on the compact space. This 
preserved the SU(2) x SU(2) isometry group of T 1,1 . There are various assumptions one could 
make to relax this condition. For instance, the probe discussion in [7J indicates that the antibranes 
embiggen to wrap an an S 2 C S 3 in the T 1 ' 1 . In this circumstance, the global symmetry is reduced 
to SO (3). Finding solutions with this reduced symmetry may be tractable. 

• Most ambitiously, one would like to find solutions for D3-branes that are valid over the whole 
KS throat. As a first step, one could find another class of perturbative solutions valid in the "near 
tip" region (where the leading order metric and background fields are very simple), and then match 
them at some intermediate r to the large r solutions we have presented in this paper. It is only by 
actually finding the solutions that extend to the IR and are built on "real" brane sources, that one 
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will be able to determine the correct quantization of the parameter S that appears in the linearized 
solutions. Similarly, only by extending our solutions to the tip is it likely that we can see the 
perturbative instability found in [7] , that sets in if one leaves the regime p <C M. 
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Appendix A: An expression for total mass 

For spacetimes asymptoting to flat space, one may define the ADM mass using the falloff of the 
metric at infinity; a similar quantity exists for perturbations of AdS space. For the Klebanov- 
Strassler cascade solution which asymptotes to neither, one needs a more general definition of the 
mass of a spacetime. We make use of a quantity formulated in more general spacetimes by Hawking 
and Horowitz [22], which reduces to the familiar cases for flat space and AdS space. Another useful 
reference is [34 J. 

For the total mass of a spacetime with a timelike Killing vector, we have 

E = --L- f N{K - /Co) , (81) 

where the norm of the timelike Killing vector is —Af 2 , K, is the trace of the extrinsic curvature and 
/Co is the extrinsic curvature of the zero-mass spacetime with the same asymptotic behavior and 
Gd is the D-dimensional Newton constant. The integral is taken at a fixed time, over a surface at 
spatial infinity; the extrinsic curvature represents the embedding of this surface in the fixed-time 
space. For a D-dimensional spacetime metric of the form 

ds 2 = -N(r)dt 2 + f{r)dr 2 + ds 2 D _ 2 {r, Q) , (82) 

where the D — 2-dimensional metric can depend on other variables Q, we may calculate the integral 
of /C as follows. The D — 1-dimensional space at constant t is foliated by constant-r slices; the 
normalized vector orthogonal to these surfaces is called n M and is 

n r = f~ l/2 , n M = , M ± r , (83) 

and one may check that this vector field generates a set of geodesies. In these circumstances, the 
(trace of the) extrinsic curvature can be written /C = Va/^ m , and we then can show 

J d D - 2 Sl^g^iK = J d D - 2 Qn M d M (^- 2 ) = n M d M A, (84) 
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with A the D — 2-dimensional area, 

A = Jd D - 2 n^- 2 . (85) 
To calculate E, the radius r of the slice must be taken to infinity. 

Appendix B: Flat space non-BPS brane solutions 

We would like to be able to provide more evidence that the mode ([20]) does indeed correspond 
in some way to a brane/antibrane perturbation. We can do this by matching to a set of known 
solutions for a stack of branes and antibranes in flat space, and taking an appropriate near-horizon 
limit. In the process we will gain a better understanding of the vanishing of the ADM mass in AdS 
space. 

Brax, Mandal and Oz constructed a set of solutions for non-BPS, Poincare-invariant branes 
embedded in flat space in various dimensions. For the 3-brane case, a class of solutions is 

e ~ 4A = cosh(kh(r))-c 2 sinh(kh(r)), e 4B = e" 4A /_/+ , (86) 
a = -7/( C 2-i) 1/2 e 4A sinh(£;/i(r)), (87) 



where k = \Jh/2 and 

/± = 1±(£) , h(r)=)og(£) . (88) 

Here we substituted f for BMO's r to avoid confusion with the quantity we use in section [3j These 
solutions correspond to a combination of D3-branes and D3-branes sitting at r = but asymptotic 
to flat space. We have set to zero a parameter c\ in [T9] which corresponds to a decoupled dilaton 
mode; since the D3- and D3-branes don't source the dilaton we chose to set its independent modes 
to zero. 

Calculating the ADM mass and integral of F§ and normalizing them to the BPS case of N 
D3-branes, 

e -4A = e 4B = a -l = 1 + R_ ? £4 = Angs ^ ^ (g9) 

we find we can identify the BMO solutions ( 1861) as corresponding to N+p D3-branes and p D3-branes 
with 



2r]r 4 k^Jc 2 2 - 1 = iT = 4vr^iV , 2f A kc 2 = 4>irg 8 (N + 2p) . (90) 

Thus the parameters of the solution r and c 2 can be traded for the total D3 charge and total mass; 
we put tildes on N, p because they do not yet coincide with the parameters N, p from our AdS 
solutions. We see that the BPS limit of p — > 0, N fixed is 



c 2 -> oo , r 4 -> , 2kf A ^c 2 2 - 1 = R A fixed , (91) 
and the expansion in p/N is equivalent to an expansion in l/4c|: 



l = C2 -V C '- 1 = ^ + ^ + (92) 
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Hence the solution can be rewritten in terms of the scale R 4 = 4irg s N associated to the background 
N D3-branes and the small parameter p/N associated to the p perturbing D3/D3-pairs: 



= l + (l + *>* + ** + l>£ + £ ><&), (93, 



iV'r 4 ' N r 8 5Nr 12 ' ~ y N 2 ' 
4B , , h , 2p^ 4 , 8p^ 8 , 2p R 12 ( 



V N r 4 5 N r s 5 N r 12 K N 2 

This is an expansion in p/N, though it would break down for r/R small with p/N fixed. The 1/r 8 
term seems to be a 1/r 4 correction to the R 4 /r 4 warp factor. This would correspond precisely to the 
dimension-four operator vev like (T^) which we did not see in the solution of section [2j However, 
these solutions for D3/D3-pairs perturbing a stack of D3-branes do not precisely correspond to the 
solution we are interested in, since they asymptote to flat space instead of AdS space. 

We can take an appropriate limit of the original solutions (1861) to achieve this, however, essentially 
corresponding to letting c 2 go to infinity such that the cosh term is dropped. Formally we can achieve 
this by rescaling the coordinates and r such that e~ 4A and e 4B acquire overall factors: 

e~ 4A = Z[cosh{k h(r)) - c 2 sinh(A; h(r))} , e 4B = e~ 4A f^f + , (95) 

and then scaling the constant Z — > with Zc2 fixed. This has the form of a perturbed AdS solution 
where the solution (l9~3l becomes 

e -« = ^ + !JL^ + 0( £1 ) (96) 
r* 5 N [ N*> ' 1 ' 

r* 5 N [ N*> ' 1 ' 

Here we defined R 4 = ZR 4 and p/N = (p/N)/Z 2 with R 4 and p/N constant as Z —>■ 0. We notice 
immediately that the 1/r 8 correction disappears along with the overall "1" in the near-horizon limit; 
there is no way to take the limit so as to preserve this correction. Moreover, it is straightforward 
to show that the perturbations in (IDIjI) are of the form 

5B = B r~ 8 , 5 A = -25B , (98) 

which is precisely the perturbation we calculated in section [2J with 

R 8 v 

Thus we see that the near-horizon limit of the solution around a stack of D3-branes perturbed by 
a small number of D3/D3 pairs is precisely the (D$ mode we already identified, with the coefficient 
(I99p . as long as a suitable scaling limit is taken. This confirms the interpretation of the Og solution 
as being (a limit of) a set of perturbing antibranes all the way down the AdS throat. Moreover 
the scaling limit to get to AdS sends the number of antibranes present in the flat space solution to 
zero, explaining the vanishing of the total mass of the AdS solution. 



Appendix C: Ricci Tensor 

In this appendix we collect expressions for the components of the Ricci tensor Rmn and Ricci scalar 
R for the metric (13T|) : 

R rr = -^-8a' 2 ---b' 2 - — -4c' 2 + a'(-— + b f + 4c')+a"-b''-4c n , (100) 
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R„u = -ri^r 2 e ia (a + rb' + Arc' + ra'(5 + rb' + Arc') + r 2 a") , (101) 
R^a = e 2a g^ a (Ae 2b - 4c -r 2 b' 2 -rb\5 + Arc') + ra'(5 + rb' + Arc') + r 2 a"-r 2 b"), (102) 
R M2 = e 2a g^ H (Ae 2b - 4c -r 2 b ,2 -rb , {5 + Arc , )+ra , {b + rb' + Arc , )+r 2 a''-r 2 b"), (103) 

R^ = ^(8e 4fe cos 2 (^) + 18e 4c sin 2 (^) - 6e 2b+2c sin 2 (^) - 2e 26+4c r 2 cos 2 (^)&' 2 

lo 

-15e 6c r sin 2 (^)c' - 12e 6c r 2 sin 2 ^)^ 2 + e 4c r(2e 2fe cos 2 (^) + 3e 2c sin 2 (^))a' 
(5 + rb' + 4rc') - e 4c r6'(10e 26 cos 2 (^) + r(8e 26 cos 2 (^) + 3e 2c sin 2 (^))c) 
+2e 2fe+4c r 2 cos 2 (^)a" + 3e 6c r 2 sm 2 (^)a" - 2e 2b+4c r 2 cos 2 (^)6" 
-3e 6c r 2 sin 2 (^)c"), (104) 

-2c 

R g . g . = (_2e 2^, + 6e 2c -e 4c r(5 + r6 , )c / -4e 4c r 2 c ,2 + e 4c ra'(5 + r6 , + 4rc') 

6 

+e 4c rV - e 4c r 2 c") , (105) 
i? = -2e 2a " 4c (2e 2b - 12e 2c + 10e 4c + 4e 4c rV 2 + e ic r 2 b' 2 + 20e 4c rc' + 10e 4c rV 2 

+e 4c r6'(5 + Arc') - e 4c ra'(-3 + rb' + Arc') - e 4c r 2 a" + e Ac r 2 b" + Ae Ac r 2 c"). (106) 

where a = ip, <f>i, % = 1, 2 and primes denote derivatives with respect to r. 

Appendix D: Matching to ABY Solutions 

Linearized solutions to the fluctuation equations around the KT background for an equivalent set 
of fields were considered by Aharony, Buchel and Yarom [IB] . They consider a considerably more 
general ansatz for these fields; one can show that our solution fits inside their framework. 

We use the Einstein frame of [17] . gf^ = S , y 2 e _ *' 2 ^?" ng , while the authors of pj5] use g^„ Y = 
e~* ' 2 <7u t " ng , so the metric we use here is related to theirs by g^ v = gl^ 2 g^ Y ■ The radial variable r 
is given in terms of their p by r = 1/p. The various functions of the solution are related as follows, 
with our notation on the left: 

e" 4a <- g s h, e 2b <- / 2 , e 2c <- / 3 , A; <-> fc, (107) 

while their 4D metric function Gij becomes Grjij for some function G(r) due to Poincare invariance, 
and then we have chosen to use coordinate invariance to set G — 1. Moreover we have the parameter 
relations 

9s <-> Po, M ^ P, iV <-> AT . (108) 

The ABY equations of motion and linearized solutions to order p 4 (1/r 4 ) are presented in their 
Appendix A. They look formidably complicated, but simplify enormously upon imposition of four- 
dimensional Poincare-invariance. Additionally, many of the parameters of their ansatz (including 
all multiple powers of logs) can be removed using a gauge redefinition 



P^ P 



1 + P 2 (^20 + $21 logp) + P 4 (^40 + <*41 logp + #42 log 2 p + S 43 log 3 p , (109) 



where the 5s are parameters of the transformation. Six parameters can thus be gauged away 
explicitly, and then a large number of other parameters are set to zero by the equations of motion. 
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Our solution precisely coincides with theirs with the two parameters S, (p corresponding to their 
modes a 4 ' and p 4 ' : 

S <-> a 4 ' , <-> p 4 ' , (110) 
and with the other nonzero coefficients determined in terms of these as 

p (4,l) = 3a (4,0) ) (nl) 

f^ 4 ' 1 ) = -ip 2 p a( 4 '°) (112) 

= ±(2K + l3P 2 p )a^-±P 2 Po p^ (113) 
64 lo 

#(4,1) = _ 3 p2 poa (4,0) (n4) 

= -(6K + llP 2 Po )a^--P 2 PoP ^°K (115) 
8 2 
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